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Abstract 

A calculation of the bulk viscosity for the massive Gross-Neveu model at zero fermion chemical 
potential is presented in the large- N limit. This model resembles QCD in many important aspects: 
it is asymptotically free, has a dynamically generated mass gap, and for zero bare fermion mass it is 
scale invariant at the classical level (broken through the trace anomaly at the quantum level) . For 
our purposes, the introduction of a bare fermion mass is necessary to break the integrability of the 
model, and thus to be able to study momentum transport. The main motivation is, by decreasing 
the bare mass, to analyze whether there is a correlation between the maximum in the trace anomaly 
and a possible maximum in the bulk viscosity, as recently conjectured. After numerical analysis, 
I find that there is no direct correlation between these two quantities: the bulk viscosity of the 
model is a monotonously decreasing function of the temperature. I also comment on the sum rule 
for the spectral density in the bulk channel, as well as on implications of this analysis for other 
systems. 
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I. INTRODUCTION 



Transport coefficients are essential inputs to describe the space-time evolution of systems 
not far from equilibrium. During the last years there has been a very active effort to analyze 
them from both the theoretical and phenomenological points of view in the context of heavy- 
ion collisions, condensed matter physics, astrophysics and cosmology. The calculation of 
transport coefficients in quantum field theory at intermediate and strong coupling is still a 
challenge from both the analytical and the numerical points of view. Due to their intrinsic 
non-perturbative nature, even in weakly interacting theories a resummation of an infinite 
number of diagrams is needed in order to obtain the leading-order result. In the strongly 
coupled regime, the most prominent method available is the AdS/CFT correspondence, 
although it is only applicable to a limited class of field theories. On the other hand, lattice 
simulations are still not accurate enough regarding the calculation of spectral densities, and 
the introduction of a finite quark chemical potential makes things even more difficult because 
of the sign problem. 

It was recently conjectured, based on a sum rule for the spectral density of the trace 
of the energy-momentum tensor in Yang-Mills theory [I], that a maximum of the trace 
anomaly near the critical temperature might drive a maximum for the bulk viscosity near 
that temperature. The corresponding sum rule was later corrected in [2], and the ansatz 
for the spectral density used to extract the bulk viscosity questioned [2-4]. Since the trace 
anomaly measures the breaking of scale invariance in a system, and the bulk viscosity £ 
essentially represents the difficulty for a system to relax back to equilibrium after a scale 
transformation, it seems in principle reasonable to think that ( would be maximum when 
the breaking of scale invariance is maximum. 

In heavy-ion phenomenology, bulk viscosity has usually been neglected because it is ex- 
pected to be much smaller than the shear viscosity even at temperatures not very high [5]. 
However, as suggested by the analysis of [1], non-perturbative phenomena responsible for 
the main contribution to the trace anomaly near T c could also produce a significant increase 
in the bulk viscosity. In this paper I will present an explicit calculation in the massive Gross- 
Neveu model in 1 + 1 dimensions, where the correlation between trace anomaly and bulk 
viscosity can be accurately tested. I will not try to give an estimation for the absolute value 
of ( in QCD near the phase transition though; as we will see this model is not suitable for 
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that purpose. There are several works analyzing this issue employing different approaches 
(see for instance [5—12] and references therein), but still the order of magnitude of the bulk 
viscosity near the crossover temperature is uncertain. 

In 1+1 dimensions, transverse flow of momentum is not possible, and the bulk viscosity is 
the only viscous coefficient present to linear order in gradients. In this paper I will analyze 
only finite temperature effects, considering a vanishing fermion chemical potential, thus the 
thermal conductivity will be zero in this case. Therefore, the only constitutive equation 
relevant for us is 1 

with (T 11 ) the non-equilibrium expectation value for the spatial component of the energy- 
momentum tensor, P eq is the pressure in equilibrium, and u 1 is the fluid velocity. The bulk 
viscosity can be in principle calculated perturbatively in field theory [13]: 

C oc lim , (2 ) 

where pbuik is the spectral density corresponding to the thermal propagator (Tff(t,x)T%(0)). 
Here though, I will use a kinetic theory approach, which should be equivalent to the dia- 
grammatic one in the perturbative (and dilute) regime [14]. 

This paper is organized as follows. In Sections II and III, I review well known properties 
of the massive Gross-Neveu model at zero and finite temperature, and I prove the breaking 
of integrability in the large-iV limit when a mass term for the fermion field is explicitly 
introduced. Then in Section IV, the calculation of the bulk viscosity within kinetic theory is 
presented. In Section V, I comment on sum rules and implications of the previous analysis 
for other systems. Finally in VI I summarize the main conclusions. There is also the 
Appendix A, where the result of factorization for fermion loops in 1 + 1 dimensions is 
derived, and Appendix B where the reader can find some details on the calculation of the 
inelastic scattering amplitude. 

1 I use the metric g = diag(+l, — 1). 
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II. VACUUM PROPERTIES OF THE MASSIVE GROSS-NEVEU MODEL 



Let's consider the Gross-Neveu model [15] with an explicit bare mass for the fermion 
field: 

N 2 / N \ 2 

C = J2^i^a + j [^Ma -Nm\ . (3) 

a=l \o=l / 

Since we are interested in studying the large N limit of the model, in order for the per- 
turbative expansion in powers of 1/N to be sensible, the bare coupling constant must be 
re-scaled, g 2 = X/N, with A being constant as N — > oo. Also, it is convenient to introduce 
an auxiliary field a to properly classify the different Feynman diagrams according to their 
topologies and power counting in 1/N [16]: 

N 1 N 

C = J ^2^ a i^ a - -a 2 - ga^2i} a ^p a + Nmga . (4) 

a=l a=l 

Clearly, the introduction of this field does not affect the dynamics of the system because 
its equation of motion is simply a = Ngm — g tp a ipa- In terms of the auxiliary field, the 
discrete chiral symmetry then corresponds to the simultaneous transformations ip ^ 1^ 
and a i— > —a. 

In 1 + 1 space-time dimensions and in the large- N limit, this model shares many impor- 
tant features with massless QCD in 3 + 1 dimensions: it is renormalizable, asymptotically 
free, classically scale invariant (for zero bare fermion mass), it has a dynamically generated 
mass gap which manifests as a peak in the trace anomaly, and in vacuum undergoes an 
spontaneous breaking of the discrete "chiral" symmetry 2 ip h> 75^. 

As we will see in the next subsection, the introduction of this bare mass m is a simple way 
of allowing the system to relax back to thermodynamic equilibrium after a small perturbation 
in the distribution of momenta. In addition, the bare mass also suppresses the density of 
kink-anti-kink configurations in the thermodynamic limit and makes the mean-field 1/N 
expansion well defined [17, 18]. 

To leading order in the large- N limit, only one counter-term is necessary to renormalize 
all the divergences, SC = 8m a a 2 /2, which essentially amounts to a renormalization of the 
coupling constant. The effective potential for the classical field a c is obtained using standard 

2 In massless QCD instead, it is the continuous chiral symmetry SU(iVf )a x SU(7Vf )y which is spontaneously 
broken in vacuum down to SU(-/Vf)v- 
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techniques and renormalized imposing the condition d 2 V e s(a c ) / dcr 2 \ Uc=ao 
the renormalization scale. This fixes the counter-term to be 

9 2 N 



1 [15], with do 



2tt 



(5) 



with A an ultraviolet cutoff. Then, the leading-order renormalized effective potential is 



Kff(^c) 



1 2 at 9 2 Na 2 c 

2 - + ^5 



In'* 



(6) 



This is a Mexican-hat potential (tilted by the mass m) with a non-zero mass gap Mq deter- 
mined by the condition 



da. 



M - g 2 Nm + In — - 2 = . (7) 

a c =M 0/9 27r [ \9°oJ 

If we define <p c = gcr c , we can now use (7) to write the effective potential in a scale- 
independent form: 

N<f>; 



Kl(0c) = Nm(f) c 



2Mr 



-1 + 



47T 



hi 



M n 



- 1 



(8) 



As shown in Fig. 1, for m = the discrete chiral symmetry is spontaneously broken by 
choosing as vacuum one of the two minima. For small enough values of m, the potential still 
has two minima, whereas for larger m one disappears and the other one becomes deeper (I 
use the units M = 1). 
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FIG. 1: Effective potential of the classical field 4> c for different values of m. 



The effective potential (6) satisfies the renormalization-group equation 



0-q 



V^{a c ,g,a ) = 



(9) 
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which implies 



g 3 N/2 



7T 



(10) 



1 + g 2 N/2n ' 

i.e., the theory is asymptotically free. Although the running coupling constant becomes 
arbitrarily large at low energies, the interaction between fermions is also suppressed by 
powers of 1/N, thus in the large- N limit we are still able to probe the low-energy regime of 
the theory. 

The leading-order contribution to the self-energy of the a field corresponds to the diagram 
depicted in Fig. 2. In Euclidean space, the expression for the (renormalized) a propagator 



m vacuum is 



1 n y M 2tt 



f3{P 2 ) In 



f3(P 2 ) + 1 



(11) 



with f3(P 2 ) = a/1 + 4Mq / P 2 a phase-space factor, and P = (pi,P2)- 




FIG. 2: Self-energy of the a field to leading order, O(N ). 



In the next subsection, I show how the first term in (11) breaks the integrability of the 
model in the large- N limit. 

A. Breaking of integrability 

The Gross-Neveu model (without the bare mass) is an integrable quantum field theory 
[19, 20], this implies the existence of an infinite number of conserved charges and 1 + 
1 dimensions the factorization of the S-matrix in terms of binary collisions, so inelastic 
processes have vanishing scattering amplitude. Since in 1 + 1 dimensions binary collisions 
cannot modify the distribution of momenta, integrability then prevents momentum transport 
in this system. Consequently, the bulk viscosity of the Gross-Neveu model is infinite. After 
including the bare mass in the model, this factorization in terms of binary collisions no 
longer happens, and hence it renders the bulk viscosity finite. 
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FIG. 3: Leading-order contribution in the large- A?" limit to the inelastic process 12 — > 3456. 



To see this, consider the leading-order diagrams corresponding to the inelastic process 
2 —7- 4 in Fig. 3. As it was shown in [20], the fermion loop of Fig. 3(g) factorizes into tree 
diagrams corresponding to all the possible ways of cutting it (Figs. 4 and 5). One particular 
cut is depicted in Fig. 5. From the result (A21) derived in Appendix A, it is easy to see 
that the four-point amplitude and the factor —T in Fig. 5 cancel out giving a —1 factor. 
Hence, the diagram of Fig. 5 exactly cancels (when m = 0) the one of Fig. 3(a), and the 
same for the rest of diagrams. If we now introduce the mass m, from (11) we see that this 
cancellation cannot happen, so the total inelastic amplitude is now oc itl/Mq to leading order 
in 1/N. This proves the non-integrability of the massive Gross- Neveu model in the large- N 
limit. 
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FIG. 4: Diagram of Fig. 3(g) expressed in terms of different cuts according to Eq. (A19). 
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FIG. 5: Factorization in terms of tree diagrams. 



III. THE MODEL AT FINITE TEMPERATURE 



The thermodynamic properties of this model have been studied in detail in many papers, 
see for instance [17, 21-24] and references therein. In this section, I am simply going to 
review leading order-results in the mean-field approximation, which are relevant for the 
later analysis of the bulk viscosity. 

The leading-order renormalized effective potential at finite temperature is 



Vj(<r c ;T) = -a 2 c -gNma c + 



Air 



In 



0-Q 



The thermal mass gap is defined by 



der. 



m 
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dk In ( 1 + e 
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dk 



n F (E k 
Eh 



(12) 







(13) 



where (7) has been used, E k = ^Jk 2 + M(T) 2 , and n-p(x) = (exp(x/T) + 1) 1 is the Fermi- 
Dirac distribution function. In Fig. 6, I plot the fermion mass gap as a function of the 
temperature for different values of m. For the case m = 0, the mass gap vanishes at the 
temperature T c ~ 0.57Mo, indicating restoration of the discrete chiral symmetry. This 
is however an artifact of the mean-field approximation; the chiral symmetry is actually 
immediately restored at T = + due to kink-anti-kink configurations 3 . Nevertheless, as 
mentioned above, the introduction of a finite bare mass suppresses these kink-anti-kink 



3 This restoration must happen in order to be consistent with the Mermin- Wagner theorem [25]. Note 
however that the phase transition at the indicated critical temperature occurs and is correctly reproduced 
in the mean-field approximation if the size of the system is kept finite and the limit N — ¥ 00 is taken first 
[24]. 
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configurations in the thermodynamic limit, and therefore we can approach in the mean-field 
approximation the curve m = as much as we wish provided we keep m finite 4 . 




0.0 0.5 1.0 1.5 20 2.5 3.0 

T/M 



FIG. 6: Thermal mass gap of the fermion field as a function of the temperature for different values 
of the quotient rajM^. 

The pressure is immediately obtained from the effective potential: 



P 



,.,«,T)/ S ;T)^Mm(l-|!)-« 

oo 

J dk ln(l + e-^W/ r ) 



V M J 



2NT 



(14) 



and the "bag pressure" is 



. iVM? ( 1 m , 



(15) 



Entropy, energy density, specific heat, speed of sound, and trace anomaly are calculated 
from the pressure using the thermodynamic relations 



A 



dP 
df 

e-P + 2P h 



e = Ts-P = T 
T 



d (P 



dT \T 
d [ P-P h 
dT 



de ds 2 dP 

gf= Qf ' c s = 



(16) 



T 2 dT \ T 2 

These are plotted in Figs. 7, 8, 9, 10, 11, 12. We see that the trace anomaly has a very 
pronounced peak right at T c for m = + , which will allow us to study the possible correlation 

4 Here it is important to emphasize that in our calculations the large-A^ and thermodynamic limits are 
taken first keeping m finite, and afterwards we study the limit m — >• + . 
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with the bulk viscosity. For m = + , above T c the pressure corresponds to an ideal gas of 
massless fermions: 



tcNT 2 ttNT 2 A 2P h 

P = — ^— , e = P , s = , c v = s , c s = l, A= — 



with P b = iVM 7(47r) 




(17) 



FIG. 7: Pressure as a function of the temperature for different values of ui/Mq. The color code is 
the same as in Fig. 6. 
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FIG. 8: Entropy density. 



In order to calculate dynamical quantities, it is convenient to shift a h-> M(T)/g + a, so 
tadpole diagrams vanish and have not to be taken into account. The a propagator to leading 
order and at finite temperature, calculated from the diagram of Fig. 2 in the Imaginary-Time 
Formalism and continued to real frequencies, is 



10 




11 



e - P + 2P h 




FIG. 12: Trace anomaly. 



— [Z\ ff (ic n ^ u + iO\p)}- 1 = - + - I + 9{s - 4M 2 )]/3( S ) In 



/3(s) + 1 



0{s) - 1 



+9(s)9(4M 2 - s)2B(s) arctan f^yj } + PV J 



1 \\ ™r f dk n F (E k ) s/3(s) 2 (2pk- s 



7i E k {2pk -s) 2 - AE 2 k u 2 



- i sgn(a;)/3( S ) - 4M 2 )[1 - n F (e + ) - n F (e_)] + 0(-s)|n F (e+) - n F (e_)|} , (18) 

where M = M(T), /3(s) = ^1 - 4M 2 /s, S(s) = ^MVs - 1, and e± = |w ± p/3(s)|/2. 

The first term in (18) is responsible of breaking the integrability of the model also at finite 
temperature, which follows from the result (A19) in Appendix A in the same way as for the 
vacuum case analyzed in the previous section. We realize that the breaking of integrability 
is now controlled by the factor m/M(T), instead of ttl/Mq. Interestingly, for T > T c , the 
limit of m/M(T) as m — > + is not zero, but a (temperature-dependent) constant. Thus, 
the scattering amplitude for inelastic processes is always finite when m = + for T > T c . 



IV. KINETIC THEORY APPROACH 



The massive Gross-Neveu model is a non-confining theory and, as we have seen in the 
previous sections, the interaction between the fundamental fermions is suppressed by powers 
ofl/N, hence in principle it seems reasonable to adopt a kinetic theory treatment to analyze 
the transport properties of this system in the large- N limit. Alternatively, one could formally 
work out the resummation of an infinite series of ladder and chain diagrams contributing the 
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spectral density of the energy-momentum tensor when the external frequency goes to zero. 
However, it is known that this resummation leads to solving an integral equation which 
coincides with the Boltzmann equation in the effective kinetic theory describing thermal 
excitations in the system [13, 14, 26-28]. 

I consider that the kinetic theory approach is simpler, and I will employ it for the cal- 
culation of the bulk viscosity in this paper 5 . I am going to follow essentially the previous 
works [5, 13, 14, 27, 30] so, although I try to keep the discussion self-contained, the reader 
is referred to these papers for additional details. 

In order to obtain the bulk viscosity, we need to determine the statistical average of 
the energy-momentum tensor of the system in a cell of fluid for a small departure from 
equilibrium. In kinetic theory, this average is [31] 

oo 

T^(t, x) = £ J ]flff A (t, x, k) , (19) 



-oo 



where f A = f A (t,x,k) is the non-equilibrium distribution function, A is a collective index 
denoting the fermionic or anti-fermionic character and the flavor of the corresponding particle 



species, Ef- = a/ M(T) 2 + k 2 , and k = (Ek, k) is the canonical momentum (the underline 
emphasizes that it is on-shell). 

The Boltzmann-Uehling-Uhlenbeck equation determines the space-time evolution of dis- 
tribution functions for dilute systems due to the change in the number of particles of type 
A produced by collisions in the fluid. In 1+1 dimensions it reads 



d + k_d\ f A Of 



dt E k dx J dt 



A 



df 



A 



dt 

gam 



1 

loss 



CM , (20) 



with / = ({f A }) a column vector containing the distribution functions for every type of 
particle. Considering only the leading-order elastic and inelastic processes in the large- N 
expansion, which in our case are 12-B-34 or 1 2 3 <H- 4 5 6, and 1 2 <H- 3 4 5 6 or 1234-H-56 

5 Strictly, due to infrared divergencies characteristic of low-dimensional systems, this calculation is valid 
in 1 + 1 dimensions only in the limit N — > oo, where the long-time tail in the energy- momentum tensor 
correlator ~ i" 1 / 2 becomes negligible [29]. Otherwise, for N finite, the bulk viscosity of the massive 
Gross-Neveu model would be infinite. 
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respectively, the collision term is given by 



cm 



dpi 



E 



2nE 1 ^ J 2tt 2tt 

B > C > D 2,3,4 



dpi dp 2 



i t-iA,B;C,D j qA,B;C,D i rC fDf-\ MVl ^ 
Qi 12->34 ^ D (l)2^34 L/3 J A \ l ~ J 1 !\ l ~ J 2 ) 



B,C,D,E,F Lr 



dp 5 dp 6 



dr t 



E,F;A,B,C,D qE,F;A,B,C,D 



2tT 27T ^56^1234 



56-s>(l)234 

2,3,4,5,6 

ErFn fCui /-D\ f AfB tC fDf-\ tE\(-\ rF\ 



x (i - /^a - /^a - / 3 G )(i - m - f?f?ftf?(i - fno- - fi)] 

dpi dp 2 



+ 



2tt 2tt 



Ar A,B;C,D,E,F T Q AB;C,D,E,F r rC rD rE rF(i fAWi /-Si 
ui 12^3456 ^ ° (1)2-5-3456 W3 J4 /5 I 1 JlJl 1 J2 ; 



2,3,4,5,1 



tf/f (1 - /?)(! - / 4 D )(1 - /?)(! - / 6 F )] + / 



dpi dj?2 d^ 3 r A,B,C;D,E,F r 2 Q A,B,C;D,E,F 
27T 27T 27T 123 ^ 456 ^(D23^456 



2,3,4,5,6 

x [/f / 5 B /6 F (i - - / 2 B )(i - / 3 C ) - tffifSQ - /f )(i - / 5 E )(i - h F )] 



(21) 



where the sum over indices runs over all the possible configurations of fermion, anti-fermion, 
and flavor states. The symmetry factors (to be specified later) £(1)24+34, 5'(i)2«>3456) 5W-Ki)234, 
and 5 , (i)23<- 5 .456 avoid the double counting from relabeling of momenta for identical particles 
(except for the particle denoted as '1') in the integral and considering equivalent processes 
after summing over all the fermion types. The transition rate for an arbitrary process a — > (3 
is given in terms of the scattering amplitude M. by [32, 33] 

dp/s 



dr{a —t/3) = L 



l—Nci 



H(2K)- 1 



n 



(2tt)2% 



|A%^/?)| 2 (27t) 2 <* (2) (E^-E 



Pf>, 



(22) 



where L is the size of the system (although we consider the limit L — > oo), and N a the 
number of particles in the initial state. 

In order to obtain an expression for the bulk viscosity we need to solve (20) for small 
departures from equilibrium. To do it, we first write 

f A (t, x, k) = f*(t, x, k) + 5f A (t, x, k) , (23) 

where 5f is small, and the fermion or anti-fermion distribution function at equilibrium for 
zero chemical potential is 

1 



icq (A x i hy 



Q 0k-u + 1 



(24) 
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with = T(t,x) the local temperature, and u^(t,x) the velocity of the corresponding 
fluid cell. Expanding the left-hand side of (20) in the local rest frame (it 1 |i. r .f. = 0) to linear 
order in spatial derivatives, we obtain 6 

du 1 



dt E k dx 1 



l.r.f. 



-n F {E k )[l - n F {E k )](3 



M dM\ 2 k 2 

Eh — i c: 

E k dT ) Et 



dx 



(25) 



where c s is the speed of sound in the fluid. 

Consequently, the deviation from equilibrium can be written in the form 

du 1 



(26) 



StfLi. = -Pn v (E k )[l-n v {E k )]Bi 

with = B A (\k\) some dimensionless function to be determined by solving the integral 
equation obtained after the previous linearization of both sides of (20): 



pI-cUeI-tm 



AM 
dT 



2(1 -n F 



E 

B,C,D 



n 

i=2 



dpi 



(27r)2Ei 



Sfn 2 ^u(^) 2si2) (Pi +P2-P3- PA)n F , 2 (l - n F)3 )(l - n FA )(Bf + B* - B% - Bf) 



'(l)2->-34 

+ E 

B,C,D,E,F 



n 

i=2 



dpi 



(2ir)2E i 



I . A A,B,C,D I Win 
\M E y (P5,P6;Pl,P2,P3,P4)| S, 



56-> (1)234 



x (27r) 2 (5 (2) (p5 + p 6 - pi - p 2 - P3 - P4)n F ,2nF,3«F,4(l - n F}5 )(l - n F , 6 ) 



x (B? + B* + B% + B? - Bf - Bf) + |A^' F (Pi,P2;P3,P4,P5,^ 
x 5 {2 \pi +P2-P3-P_4-P5 -pp)n Ft2 (l - ra F)3 )(l - n F)4 )(l - n Fj5 )(l - n F)6 ) 
x (Bf + Bf - S 3 C - B° - Bf - Bf) + {M^^^p^^S^^ 
x (27r) 2 <5 (2) (pi +p_2+P3-Pi-pp - Pp)n Fj2 n F ^{l ~ n Fj4 )(l - ri Fj5 )(l - n F , 6 ) 
x + B b + go _ b d _ b e _ bF) 



1 2 aAB-.CDEF, 



(27) 



We can interpret the right-hand side of (27) as the action of a linear operator c <o over a 
function in the space of solutions of the transport equation, and we split this operator into 
two terms, ^ = corresponding to elastic and inelastic processes respectively. At this 



6 Here we make use of the thermodynamic relations dT/T = dP/(e + P), = dP/de. Also, from the 
conservation of the energy-momentum tensor, d fl T llu = 0, applied to leading order to the perfect fluid, 
Tp f = —Pg^ v + (e + P)u^u w , we derive the relations in the local rest frame (9 ft u°|i. r .f. = 0): de/dt — 
-(e + P)du 1 /dx, du^/dt = -(e + P)~ 1 dP/dx. 
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point, an important simplification is in order. Since the source term in (25) is invariant under 
charge conjugation, and the theory is symmetric under 0(iV)-flavor rotations (this symmetry 
cannot be broken in 1 + 1 dimensions [34]), then the departures from equilibrium are the 
same for all the particle types, i.e., Bf = S(|&|). Furthermore, the ^-function in the elastic 
2 — y 2 term of the collision integral implies in 1 + 1 dimensions that the final set of momenta 
are the same as the initial, i.e., 2 — > 2 elastic collisions in 1 + 1 dimensions cannot relax 
back to equilibrium a perturbation in the distribution of momenta. Thus, %(2 — > 2) = 
and therefore, to leading order in the large- N expansion, = ^ e i(3 — > 3) + % n (2 -h- 4). 

Once we know from (26), (1), and (19) it is straightforward to obtain the bulk 

viscosity: 7 



oo 

dM 



C = /3 ^/ ^- k ME k )[l - n F (E k )] 



k 2 -c 2 [E 2 h -TM 

J I 'Ztt hu 

A 



B A {k) . (29) 



-oo 



The linearized version of the Boltzmann equation (27) can be written as 

\S) =<f\B), (30) 

where S(p) = p 2 -c 2 (E 2 -TM dM/dT) denotes the source term. Defining the scalar product 
of two square-integrable functions as 

oo 

(m / 7^rn F (E k )[l-n F (E k )} X A (m A (k) , (31) 

A -oo 

then the bulk viscosity is given by 

( = (S\B) = (S\tf- l \S) . (32) 

As shown in [5, 27, 30], in order to calculate numerically this expectation value for the inverse 
of the collision operator, it is optimal to do it variationally. If we define the functional 

Q\X\ = (X\S)-\(X\^\X), (33) 

7 The Landau-Lifshitz condition 

oo 

/dk 
—— (E 2 k - TMdM/dT)n F (E k )[l - n F (E k )]B A {k) , (28) 
lixE k 

— oo 

imposed to make the decomposition (23) unique, is also used here [14, 35]. 
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then the solution of (30) corresponds to a maximum in this functional, SQ/5x\ x=B = 0. 
Hence, the bulk viscosity is proportional to this maximum: 

C = 2Q max . (34) 

We now expand the solution for the Boltzmann equation in terms of a given set of n linearly 
independent functions: 

n 

= • (35) 

i=i 

Then 

n n 

Q[{h}\ = ~ 2 biCijbj ' (36) 

i=l i,j=l 

where 

Si = (0,15) , Cy = (0 4 |^|0,> . (37) 

Maximizing (36) with respect to the set of coefficients implies b = C^ 1 ^ (a tilde denotes 
matrices), and therefore 

( = g% = S'C^S . (38) 

It is important to notice, from (27), that the collision operator has one zero-mode Xe(p) = E p 
corresponding to energy conservation, i.e., ^\Xe) = 0. 8 Therefore, in order to be able to 
invert the collision matrix, it is necessary to calculate it in the vector space orthogonal to 
this zero-mode. This does not affect the result for the bulk viscosity because the source 
term is orthogonal to this zero- mode, (S\Xc) = 0. 



In addition to \c, the elastic collision operator also has the zero- mode Xn(|fc|) = 1 corresponding to the 
conservation of the total number of particles in these type of processes. However, this is not a zero-mode of 
the inelastic part of the collision integral and therefore we do not have to worry about it when calculating 
^€ . The presence of this other zero-mode, though, implies that the bulk viscosity is dominated by 
inelastic processes at very low temperatures due to Fermi-Dirac factors: 

C=(5|^- 1 |5)^ ,i (X " |tS)l l, . , forT«M . (39) 
(Xn|^i„(2o4)|x„) 

On the other hand, at temperatures close to T c , the Fermi-Dirac factors are 0(1) and both 3 — > 3 and 
2h4 processes are equally important. 
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Besides the simplifications already commented in the previous paragraphs, in order to 
obtain an explicit expression for the matrix element of the collision operator in the large- 
iV limit, it is obvious (cf. Appendix B) that the dominant scattering processes are those 
for which three different flavors participate. Moreover, the symmetry factors have to be 
specified, see Table I. 



3—^3 processes 


'S'(l)23-5>456 


£Cl£&£C y £d£ft£C 


1/12 


£a£Cl£6 £~a£a£b y £C£C£f> 


1/12 


£a£a£c £~a£a£c y £~fc£b£c 


1/12 


£a£&£?> y £C£C£d 


1/12 


£a£C£C y £fe£b£a 


1/12 


£a£&£C v £a£fe£C 


1/12 


£a£C£b v £a£C£6 


1/12 


£a£&£C v £a£fe£C 


1/12 



2f>4 processes 


'5(1)2^3456 


(1)234 


£a£fc ^_y £a£fe£C£C 


1/24 


1/12 


£a£C ^_y £a£C£()£fe 


1/24 


1/12 


£fe£C ^_y £a£a£&£C £a£a£&£C 




1/12 


£&£& ^ £a£a£C£C £a£a£C£C 




1/24, 1/12 


£C£C ^ £a£a£&£f> £a£a£b£fe 




1/24, 1/12 


£~ft£C ^ £a£a£fo£C £a£a£fe£C 




1/24, 1/12 


£C£& ^ £a£a£C£b £a£a£C£fe 




1/24, 1/12 


£a£~6 /_v £a£b£C£C 


1/48 


1/24 


£a£C /_v £a£&£C£fe 


1/48 


1/24 


£a£fe <_> £a£C£C£fe 


1/48 


1/12 


£a£C /_v £a£&£b£C 


1/48 


1/12 


£a£a £~a£a ^ £{>£&£C£c 


1/48 





TABLE I: Symmetry factors corresponding to relabeling the momenta of identical particles under 
the collision integral after summing over all particle types. Here a 7^ b 7^ c 7^ a denote only flavor 
(processes with three different flavors dominate). The momenta for each particular process are 
initially labeled according to the order indicated in the title of columns 2 and 3. The particle 
labeled with '(1)' is "distinguishable". Some cells are empty because the corresponding process 
has already been taken into account by another symmetry factor. The processes obtained by 
permutation of the three flavors, although omitted in the table, have the same symmetry factors 
and also have to be taken into account. 
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Finally, 



Cij ~N S 



n 



dpi 



(2vr) 2 5 (2) (pi +P2-P3-P4-P5- P6)n F ,in Ffi 



I 12—^3456 1 2 + 2 1-^^12— ^3456 1 2 



\ (2vr)2E 4 

x (1 - n F>3 )(l - n Pi4 )(l - n F)5 )(l - n F)6 ) 

x [0<(pi) + 0i(pa) ~ <t>iip?) ~ ^i{Pi) ~ <MPs) - 0i(pe)] 
x [0i(Pi) + - 0i(pa) - 0i(P4) - 0i(Ps) - 0j(Pe)] 

+ (27r) 2 5 (2) (pi +P2+P3-P4-P5 -P6)^F,l n F,2^F,3(l ~ ™F,4)(l ~ ^F,5)(l ~ «F,l 

[0i(Pl) + 0i(P2) + 0i(ps) - 0i(p4) - 0i(PB) ' 
X [0j(Pl) + 0j(P2) + 0i(Ps) - 0i(P4) - 0j(ps) - 0j(Pe)] } , 



X 



I 3 ' 

77 |A^i23->456| + 7: I -A^ 123^-456 1 

o 2 



0i(Pe)] 
(40) 



where .M 12-1.3456 and 12-5.3456 are the inelastic amplitudes of fermion-fermion and anti- 
fermion-fermion scattering respectively. The amplitude squared \M. 12^3456 1 2 °f anti-fermion- 
anti-fermion scattering, as well as its corresponding symmetry factor, are equal to the 
fermion-fermion ones by charge conjugation, and have already been included in (40). 

It is evident from (40) that the collision matrix is symmetric and positive semidefmite 
(positive definite in the space orthogonal to its only zero-mode). Note also that Si = O(N), 
and since Cy = 0(1/N) (cf. Appendix B), therefore ( = 0{N 3 ). 

A. Numerical results 



A particularly convenient set of functions, which becomes a basis when n — > 00, is [5 



(41) 



;i + l*l/<l*D)"- 8 ' ' 

with the thermal average (\k\) ~ \JMqT for T — > 0, (\k\) ~ T for T — > 00, and interpolating 
between these two behaviors for intermediate temperatures. This set of functions auto- 
matically incorporates the required asymptotic behavior for the solution of the Boltzmann 
equation in the bulk channel: ~ 1 for \k\ — > 0, and -B(|/c|) ~ k 2 for \k\ — > 00. 

In Fig. 13, I plot the numerical result of a variational computation of the bulk viscosity 
in the massive Gross-Neveu model with m = 10~ 2 Mo using n = 3 basis functions. It is not 
difficult to realize that ( increases exponentially at low temperatures, like ~ exp(2Mo/T), 
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due to the Fermi-Dirac factors present in C as well as in S (cf. footnote 8). This behavior 
is analogous to the case of ( for A(/> 4 in 3 + 1 dimensions [13]. 

The numerical error corresponding to considering only n = 3 basis functions (including 
the error from the numerical evaluation of integrals) is estimated to be of order 0.5% for 
temperatures around T c , increasing as we go down in temperatures, being ~ 60% for T = 
O.lMo, which indicates that the basis (41) is not the best choice at those temperatures. By 
considering n = 9, the precision can be improved to ~ 20% at T = 0.1M . However, due 
to the exponential growth of ( at low T and since the result shown corresponds to a lower 
bound, the qualitative behavior with temperature is not expected to change significantly. 
From the numerical result we clearly see that there is no maximum in the bulk viscosity 




i i i i i i i i i i i i i i i_ 

0.0 0.5 1.0 1.5 

T/T c 



FIG. 13: Bulk viscosity of the massive Gross-Neveu model for m = 10 2 Mo, calculated with n = 3 
basis functions. The continuous line simply joins the data points. 

near T c , it is a monotonously decreasing function of the temperature. By reducing further 
the value of m, we would eventually reconstruct (continuously) a discontinuity for £ at T c . 
For infinitesimally small m, above T c the bulk viscosity would be arbitrarily small. Going 
down in temperatures, it would increase very sharply right at T c (with an arbitrarily large 
value 9 ), and it would continue increasing exponentially at very low temperatures. This is 
shown in the plot of Fig. 14. 

9 This is perfectly fine for our purpose of testing the possible correlation between the bulk viscosity and the 
trace anomaly; nonetheless, this model is not suitable to obtain for instance an estimate of the absolute 
value of the quotient (/s for QCD. 
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FIG. 14: Bulk viscosity of the massive Gross-Neveu model for m = + , calculated with n = 3 basis 
functions. The continuous line simply joins the data points. 

V. DISCUSSION 

A. Sum rule in the bulk channel 

In the paper [2] (Section IV), the authors obtained a sum rule for the spectral density 
of the two-point function involving the trace of the energy-momentum tensor in pure Yang- 
Mills theory. The derivation is essentially based on the asymptotically-free character of 
the theory 10 . In this subsection I am going to comment on an interesting aspect of the 
massive Gross-Neveu model concerning its sum rule in the bulk channel 11 . Naively, since 
the massive Gross-Neveu model is asymptotically free, we can follow the analysis in [2] and 
convince ourselves that the version of the sum rule for this particular system is simply 

oo 

(e + P)(l-a-2(e-P) = ^|dc,^M. (42) 

o 



However, there is a recent example of an asymptotically-free model for which the sum rule in the bulk 
channel has a different form from the one derived in [2], due to the fact that conformal symmetry is not 
restored at high energies or temperatures in this model [36] . Note instead that in the massive Gross-Neveu 
model, although m explicitly breaks scale symmetry, it is eventually restored at high energies because m 
appears in the lagrangian multiplied by g 2 N . 

A more rigorous and detailed analysis of the sum rule and the spectral density is underway and will be 
published elsewhere. 




m = 



J i i i I i i i I i i i L 



.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

T/T c 
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Interestingly, we notice that this sum rule does not depend explicitly on the mass parameter 
m. Since e, P = 0(N) and c s 2 = 0(N°), the LHS of (42) is 0(N). On the other hand, the 
order of <5pbuik/w depends on the frequency. For frequencies of order u> ~ M , a diagram 
contributing to <5pbuik/w a t leading order is the one of Fig. 15, which is 0(N). Hence, the 
contribution to the integral in (42) from this region of frequencies is O(N), consistent with 
the LHS of the sum rule. 




FIG. 15: One of the diagrams contributing to 5pbuik/ w at leading-order for oj ~ Mq. The crosses 
denote insertions of the trace of the energy-momentum tensor. The fermion propagator is not 
dressed. 

If we now consider lower frequencies, u ~ then a resummation of diagrams is neces- 
sary to obtain the leading-order spectral density. This is essentially because of the presence 
of pinching singularities when the external frequency becomes of the order of the fermion 
width 7f ~ ImE = 0(1/N) (cf. Fig. 16) [13, 26, 37-39], or smaller. These singularities cor- 
respond to the product of retarded and advanced fermion propagators sharing approximately 
the same momentum: 

Sret(P)S adv (P) ~ — = O(N) . (43) 
7f 

Consequently, the set of ladder diagrams depicted in Fig. 17 all contribute at leading order, 
0(N 2 ), to the spectral density in this range of frequencies. Thus, their contribution to the 
integral is again consistent with the LHS of the sum rule. 




FIG. 16: Leading-order contribution to the fermion self-energy. 

For even lower frequencies, we know from the analysis of the Boltzmann equation in the 
previous section that 3 — > 3 and 2 -H- 4 processes will eventually dominate the spectral 
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FIG. 17: A ladder diagram with an arbitrary number of rungs contributing to <5pbulk/w a t leading 
order for frequencies u ~ 1/N. Double lines denote dressing of the fermion propagator according 
to the diagram of Fig. 16. The presence of pinching singularities is crucial in this regime of low 
frequencies. 

density. Since the bulk viscosity is of order 0(N 3 ), for frequencies close to zero 5pbuik/w = 
0(N 3 ), and therefore, in order for it to be consistent with the LHS of the sum rule, this 
region of frequencies must be w < 1/N 2 . Now, in the previous section we saw that for 
temperatures T < T c , if we reduce the value of m, the bulk viscosity can become arbitrarily 
large. However, the LHS of (42) remains finite as m — > + , this implies that the region 
where ^pbuik/w is of order 0(N 3 ) has a width ~ 1/N W with w > 2 and therefore the bulk 
viscosity does not contribute to the sum rule in this regime of temperatures for any value of 
m. 

B. Other systems 

From the analysis of the previous sections for the massive Gross-Neveu model we can 
already extract some conclusions for other similar systems. Consider for instance the non- 
linear cr-model in 1 + 1 dimensions [10]. This model also shares with massless QCD the 
features of asymptotic freedom, dynamical generation of a mass gap, and classical scale 
invariance broken by the trace anomaly. The lagrangian of the model is 



Since the O(N) symmetry cannot be broken in 1 + 1 dimensions, there is no phase transition 
in this system at finite temperature. It again is convenient to introduce an auxiliary field a 
in order to analyze diagrammatically the large- iV limit: 



This model is also integrable [20]. To leading order in the 1/N expansion, the inelastic 
diagrams have the same topology as the ones in Fig. 3. Integrability is proven in an 




a = 1, . . . ,N with the condition a a = 1/g 2 . 



(44) 




(45) 
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analogous way to the Gross-Neveu model using the factorization of scalar loops in 1 + 1 
dimensions. And it can be broken for instance by introducing a term ~ n(<j) a <j) a ) 2 in the 
lagrangian, so the correlation between the bulk viscosity and the trace anomaly can be 
studied by making k arbitrarily small. 

The thermodynamic properties of the system have been studied for instance in the works 
[41, 42]. In Fig. 18 I plot the thermal mass gap, speed of sound, and trace anomaly. There 
are no discontinuities for these quantities in the limit k — > + . We then realize that the 
qualitative behavior of the bulk viscosity as we decrease k (restoring integrability) is the one 
depicted in Fig. 19. For this model, integrability is always restored as k — > + and therefore 
bulk viscosity diverges. One could also consider the sum rule, which presumably is identical 
to the case of the Gross-Neveu model, and similarly the bulk viscosity would not contribute 
to it. 




FIG. 18: Thermal mass gap, speed of sound, and trace anomaly of the non-linear u-model in 1 + 1 
dimensions to leading order in the 1/iV expansion. (•)* denotes the finite-temperature part. 




the breaking of integrability 
decreases in this direction 



Too 



FIG. 19: Qualitative behavior for the bulk viscosity of the non-linear cr-model in 1 + 1 dimensions 
in the large- iV limit. 

Pure Yang-Mills theory in 3 + 1 dimensions in the large- N c limit is also very similar to the 
massive Gross-Neveu model regarding the bulk viscosity. For low energies and temperatures, 
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interaction between glueballs is suppressed by powers of N c , being the 3-point scattering am- 
plitude ~ 1/N C [43]. In this case the bulk viscosity is dominated by inelastic processes, due 
to the presence of a zero-mode in the collision integral corresponding to particle-number 
conservation [13, 14]. This implies an exponential growth of the bulk viscosity as the tem- 
perature decreases, ~ exp(m g /T), with m g the mass of the lightest glueball. Regarding the 
sum rule derived in [2] for this theory, at temperatures below T c , the region of frequencies 
where inelastic processes dominate becomes very narrow to be consistent with the sum rule. 
This implies for instance that extracting transport coefficients on the lattice at temperatures 
below T c becomes much more difficult as N c increases [37]. In this regime of temperatures, 
a kinetic theory approach instead is more suitable in the large- N c limit to obtain transport 
coefficients. 

On the other hand, massless Nf = 2 QCD is qualitatively different from the previous 
models. At very low temperatures, the dynamics is dominated by Goldstone bosons and 
there is no exponential growth in the bulk viscosity as the temperature decreases. In addi- 
tion, since this system undergoes a second-order phase transition in 3 + 1 dimensions, the 
bulk viscosity would diverge at the critical temperature [4, 9, 35, 44]. 

VI. CONCLUSIONS 

We have seen that the massive Gross-Neveu model is non-integrable in the large- N limit, 
which allows the study of momentum transport in this system. We found that there is 
no direct correlation between the trace anomaly and the bulk viscosity in general, i.e., a 
peak in the former does not necessarily imply a peak in the latter. 12 This was already 
obtained by S. Jeon in [13] where he analyzed the bulk viscosity of massive A0 4 theory in 
3 + 1 dimensions, but since it is not an asymptotically-free theory and the scale symmetry 
is explicitly broken in that case, it remained to analyze whether a QCD-like theory could 
be qualitatively different due to the anomaly (idea originally motivated by the paper [1]). 
The use of this simple model in 1 + 1 dimensions also avoids interference from critical 
phenomena present in higher dimensions (for instance, the bulk viscosity would diverge 
near a second-order phase transition). In addition, it is a useful model to study sum rules 

12 Of course these two quantities are not independent of each other; conformal theories have zero bulk 
viscosity. 
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and transport coefficients both at zero and finite fermion density in the large- N limit. For 
instance, regarding sum rules, after a first superficial analysis, we saw that below T c the 
bulk viscosity would not contribute to the sum rule. This implies in general that it is not 
necessarily possible to extract bulk viscosity from sum rules. Further work in these directions 
is in progress. 
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Appendix A: Factorization of fermion loops in 1 + 1 dimensions 

In 1 + 1 dimensions, the Lorentz group consists only of boots and therefore fermions have 
no spin. A two-dimensional representation of the Dirac algebra is for instance 

7° = a 1 =r 7l = i(T2= (_° i J,) => {7 M ,7 1/ } = 2^, (Al) 

with a 1 the Pauli matrices, and g = diag(+l, —1). 

The general solution of the Dirac equation (i , ~f tJ- d fl — m)ip(x) = can be written as a 
linear combination of plane waves u(p) e~ ip ' x , v(p) e ip ' x , corresponding to fermions and anti- 
fermions respectively, where p° > and p 2 = m 2 . The spinors are normalized as 

u(p)u(p) = 1m , v(p)v(p) = —2m , (A2) 

and they also verify 

u(p)u(p) = i> + m , v(p)v(p) = $ — m . (A3) 

In the rest of this Appendix, I will derive the finite-temperature version of the result 
previously obtained in [45] concerning the factorization of fermion loops in 1 + 1 dimensions. 

26 



Consider the momentum integral and Matsubara sum corresponding to the fermion loop in 
Fig. 20 at finite temperature with n > 3 (which is finite in 1 + 1 dimensions): 

oo 

n ~ ^ J 2tt (K + Q^ + M 2 ••' (K + Q n ) 2 + M 2 ' 1 ' 

Wm -oo 



with K = (—uj m , k), uj m = (2m + l)nT (m G Z), = (— i/j, <&), ^ = 2r7rT (r G Z). 



P 3 




FIG. 20: A fermion loop with n > 3 external legs corresponding to the a field. The different 
momenta satisfy Qi = Qi-\ + Pi, with Qq = Q n = 0, and Ya=1 ^ = 



We first perform the Matsubara sum using the result [39] 

oo 

T F ^m) = ReS ( F ' ~ £ / 2^ ^ DiSC(F ' CUt) ' (A5) 



poles 2i 

13 



The integrand in (A4) has poles at iu m = — iz/j ± E k+f}i , and no cuts. Thus 

" = ^ S£ MsW F5^; n ^wtm- • (A6) 

-OO i_i S - =C± 3 = 1 W 

where Qji = Qj — Qi, and Kf* = (isE)~ +qi , k + qi) is on-shell, i.e., (K^) 2 = —M 2 . If we now 
make the change of variables 



13 In what follows I will be careless with Dirac indices, what makes manipulating expressions easier. However, 
we must somehow keep track of them so the hnal result has the right matrix structure in Dirac space. 
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then the integrand in (A6) becomes rational: 

— 1 oo N 



1 I f f 1 " 

L « = 27 ^ / + / <"* S ^(^))7T 

* =1 V-oo 1 / * =±1 1 



M-$ 



(0 



«)(^ 2 -l) 



j = 



(A8) 



where JE^ and Kf~ mean making the substitution k >->• k(li,qi) in and Iff^ according 

to (A7), explicitly 



Z 2 + 1 



a: 



.1/ 



w Z? - 1 



:i S (/ 2 + i),2/, 



and are the solutions of the equation 



' 



(K s (i) + Q^f + M 2 = <£> 



2Mq Jhl ±\Q 2 ji \Jl + 4MyQ 



'31 



'CO 

After partial fraction decomposition, we have 

-1 oo 



i2sM te . 2 + Q J 2 j 



5 



11 



S=±l 1=1 

n 

-EE 

s=±l i=l 

where obviously 



4 + 



.4 



.7' 



-oo 1 



-1 oo s 



3 = 1 

3 # » 



dh n F [ sm|^ )l s (/i,Q. 



-oo 1 



= lim 

i-+i%* 



And explicitly, 



.4 



1 1 

27 HO 



[M-Jji 



(y'±)' 



n 



^— (i2 S M fc>2 + - - 1^) 



(»j'±) w J*' 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 



14 The terms corresponding to the poles at U = ±1 give a zero contribution to the integral (it can be easily 
seen from the fact that the integral must finite), so we ignore them. 
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where K%- ± ^ means making the substitution /j i-> l^' s in Kf^. From Eq. (A10) it is evident 
that the momentum + Qji is also on-shell. And using the relation (A3), rotated to 

Euclidean space, we can rewrite this expression as 



4±,s 



x 



n 

*: = i 



' kiJ 



T 



±,s 



T±snQlyJl + lM*/Ql 



(A15) 



Also note that the momenta + Q k i with k ^ i,j instead, are not on-shell. Now, since 

K(ij±) and ^(y±) + Qji are both on-shell, this implies 15 



K(a+) + Qji - Kf H 



(A16) 



Then, from (A15) and (A16) we get 



ji ~ ij 



rps 

ji • 



(A17) 



Furthermore, from the property 
previous case it is not difficult to obtain 



K s=l 



, in a way analogous to the 



rpS=L rps = —l rp 

1 ji ~ ± ij = 1 j i 



(A18) 



Then, working out the integral in (All) and using the previous symmetry relations, after 
a tedious (although straightforward) simplification, one arrives to the final result 



Eiji 
At, 



i,j=l 
i # 3 



4tt 1 0(Q%)Q% 



111 



g(gl) + 1 
mi) - 1 



-2 / dk 



n F (E k ) Q% + 2k gjiil 

E k (Q% + 2kq ji , 1 y + iEl q l 2 

(A19) 



lo 



Because then (2K s {lJ±) + Q Si ) ■ Qji = 0, and (-2^ i _ 1 + Q jt ) ■ Qji = 0. Thus, 2Kf. . +) + Q jt = a(2^ + 



(y+) 



Qji) (in two dimensions there are only two linearly-independent vectors), so a — — 1 (after looking at (A9) 



and (A10)), and therefore K { 
so a = — 1, and X, 



K (ij+)-Qji- On the other hand, 2Kf. j+) +Q jt = a(-2K s Gi _^+Qji 



{ji- 



K. 



(«+) ^ ^ z 
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with /3(Q?i) = J 1 + AM 2 /Q'j i . The result of [45] is then obtained by particularizing for 
T = and performing the rotation to Minkowski space. 



We can represent the result of (A19) diagrammatically as in Fig. 21, with 




(A21) 



FIG. 21: Diagrammatic representation of the factorization of fermion loops in 1 + 1 dimensions 
into two tree graphs, as implied by the result (A19). The momenta Kuj^ \ and Kuj^ + Qji are 
on-shell. 



Appendix B: Scattering amplitudes 

As we saw in Section IV, in 1 + 1 dimensions only 3—7-3 and 2 <H- 4 scattering amplitudes 
contribute to the bulk viscosity at leading order in the 1/N expansion. Among them, in 

16 Note that the contraction of Dirac indices due to the trace in the original loop diagram of Fig. (20) is 
apparently gone, but it is actually there because the external momenta of the trees are equal at both 
sides, and therefore 

[u(q)Au(p)][u(p)Bu(q)] = [u a (q)A ab u b (p)][u c (p)B cd u d (q)} 

= u d (q)u a (q)A ab u b (p)u c (p)B cd = Ti{u(q)u(q)Au(p)u(p)B} , (A20) 
for any momenta p, q and any matrices A, B. 
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the large- N limit, processes involving three different flavors dominate because of an overall 
factor N 3 obtained after summing over all flavor types. 

By following the discussion in Section II A, we realize that the diagram of Fig. 3(g), after 
being cut in all possible ways, cancels the contribution from the diagrams 3(a-f), except for 
a term proportional to the integrability-breaking parameter, i.e., m/M(T). And this also 
happens in the case of diagrams contributing to the 3 — > 3 amplitude, which are of the same 
order as the 2 -H- 4 ones and have the same topology. 

Consequently, the fermion-fermion inelastic and three-fermion elastic scattering ampli- 
tudes are given by the sum of the diagrams of Fig. 22 and 23 respectively times the 
integrability-breaking factor with a minus sign: 



rn 



iMtt^ttft = -i - M [a) - +M (b) - +M (c) - +M (d) - + M if) 

v v H f->f fff M(T) ff ^ ffff ff->ffff ff^ffff ~ J l ff->ffff ~ J l ff^ff I 7 



LMfff-^.fff — — i 



rn 



M{T) 



M [a) +M ( ) + M (c) +M ( ) + M U) 

- /Vl fff^fff > • /Vl f f f->f f f > Jvl {f f->f ff ' JVi i ff->fff > - /Vl f f f— s-f 



f f 



(Bl) 



(B2) 



and similarly for the processes involving more anti-fermions 17 . 
For instance, the explicit contribution from the diagram 22(a) is 

LM$_# m = HP3)ST(P3+P5+P6Hpi)u(p i )u(p 2 )u(p 5 )v(p e )Dl c \p 4 -p 2 )Dl c \p 5 +p 6 ) , 

(B3) 

and analogously for the rest of diagrams. We already see from (B3), since D r ^ = 0(1/N), 
that the amplitude squared is \A4 { f _>.j f f f | 2 = 0(1/N 2 ) (and the same for the other ampli- 
tudes). 

It is not difficult to realize that the momenta in the argument of the retarded fermion 
propagators cannot be on-shell, being each p^ on-shell, therefore we can make the substitution 
Sp 1 i — y Sp(po,pi) = l/(j/>— M) in our calculations without worrying about singularities. 
After summing all the amplitudes and squaring them, we can simplify the expression using 



17 By symmetry under charge conjugation, \M { f _>.f f f f | 2 = \M.j f_>j j j f | 2 , \A4f f f_>.f t f | 2 = \M.j j j_>f j j | 2 , and 
l-^f f f-ff f f | 2 = l-^f f f->f f f | 2 - 
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FIG. 22: Diagrams contributing to the inelastic process f f — >• f f f f in the massive Gross-Neveu 
model. Here a, 6, and c denote arbitrary flavors. The diagrams corresponding to f f — > f f f f are 
obtained inverting the fermionic flow (but not the momentum) in the line which joins p\ and p%. 
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FIG. 23: Diagrams contributing to the elastic process fff— ?• f f f in the massive Gross-Neveu model. 
Here a, b, and c denote arbitrary flavors. The diagrams corresponding to fff— >■ f f f are obtained 
inverting the fermionic flow (but not the momentum) in the line which joins pi and p±. 
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(A3), the property D* ret (p) = D a ^ et {—p), and the relations for traces in Dirac space 

Tr{(jfi+M)(4+M)} = 2(M 2 +p-q) , (B4) 
Tr{(^+M)(^+M)(^+M)} = 2M[M 2 +p - (k + q) + k- q] , (B5) 
Tr{(^+M)(^+M)^+M)(/+M)} = 2[M 2 (M 2 + / • (p + fc + g) + p ■ (fc + g) + fc ■ q) 

+ (/ • p)(fc ■ q) - {I ■ k)(p .q) + (l.q)(p. k)} . (B6) 

Since the final expressions for the squared amplitudes after simplification are still very long 
and do not explicitly provide further significant information, I avoid to write them down 
here. 
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